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ABSTRACT 



We observed with the SOPHIE spectrograph (OHP, France) the transit of the HAT-P-6b exoplanet across its host star. The resulting 
stellar radial velocities display the Rossiter-McLaughlin anomaly and reveal a retrograde orbit: the planetary orbital spin and the stellar 
rotational spin point towards approximately opposite directions. A fit to the anomaly measures a sky-projected angle A = 166° ± 10° 
between these two spin axes. All seven known retrograde planets are hot jupiters with masses M p < 3 Mj up . About two thirds of the 
planets in this mass range however are prograde and aligned (A = 0°). By contrast, most of the more massive planets (M p > 4 Mj up ) are 
prograde but misaligned. Different mechanisms may therefore be responsible for planetary obliquities above and below ~ 3.5 Mj up . 

Key words. Planetary systems - Techniques: radial velocities - Stars: individual: HAT-P-6 



1. Introduction 

Spectroscopic observations during the transit of an exoplanet 
across its host star can measure the sky-projected angle between 
the spins of the planetary orbit and the stellar rotation (the obliq- 
uity) through the Rossiter- McLau ghlin (RM) effect (Holt |1893l 
Rossiter 1924; McLaughlin 1924). The occultation of a rotating 
star by a planet distorts the apparent stellar line shape by re- 
moving the profile part emitted by the hidden portion of the star. 
This induces anomalous stellar radial velocity variations during 
the transit, which constrain the sky-projected obliquity (A) and 
thus indicate whether the orbit is prograde, retrograde, or polar. 

Queloz et al. (2000) reported the first detection of the RM 
anomaly for an extrasolar planet, HD 209458b. That planet 
shows an aligned, prograde orbit, as did all of the first seven 
planets for which the RM effect was measured (as reviewed in 
Hebrard et al. 2008). These early results were interpreted as 
a validation of theories of planetary formation and evolution 
where a single giant planet migrates in a proto-planetary disk 
perpendicular to the stellar spin axis (e.g. Lin et al. 1996 ). That 
migration is expected to conserve the initial alignment between 
the angular momentums of the disk and of the planetary orbits. 

Hebrard et al. (2008) however found a first case of spin- 
orbit misalignment for the XO-3b planet, confirmed by Winn 



* Based on observations collected with the SOPHIE spectrograph 
on the 1.93-m telescope at Observatoire de Haute-Provence (CNRS), 
France, by the SOPHIE Consortium (program 10A.PNP.CONS). 
SOPHIE radial velocities are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



et al. d2009al >. Thereafter, Moutou et al. d20091 > reported a 
second case, HD 80606 (see also Pont et al. 120091 Winn et 
al. 1200 9b, Hebrard et al. 1201 01) . A dozen misaligned systems 
have now been identified, including some with retrograde or 
nearly polar orbits (e.g. Winn et al. I2009ct Narita et al. 2010a; 
Triaud et al. 120101 Simpson et al. 120101 . These unexpected re- 
sults favor alternative scenarios where close-in massive plan- 
ets have been brought in by planet-planet (or planet-star) scat- 
tering, Kozai migration, and/or tidal friction (e.g. Malmberg et 
d. 120071 Fabryc ky & Tremaine 120071 C hatterj ee et al.[2008; 
Nagasawa et al. 120081 Guillochon et al. 120101 ). Alternatively, 
it has been proposed that the orbit still reflects the orientation 
of the disk, with the stellar spin instead having moved away, 
either early-on through magnetosphere-disk interactions (Lai et 
al. 120101 ). or later through elliptical tidal instability (Cebron et 
al. 1201 11 1. Distinguishing between these mechanisms needs ad- 
ditional obliquity measurements (e.g. Morton & Johnson ROlOl l. 

Here we present spectroscopic observations of one transit of 
HAT-P-6b. This hot jupiter transits a bright F star (V = 10.5) 
every 3.8 days and was discovered by Noyes et al. (12008, N08). 
Its mass is 1.06 + 0.12 M Jup and its radius 1.33 + 0.06 R Jup . 

2. Radial velocity measurements with SOPHIE 

The August 21 th 2010 transit of HAT-P-6b was observed with 
SOPHIE. This cross-dispersed, stabilized echelle spectrograph 
is dedicated to high-precision radial velocity measurements 
(Pen-uchot et al. 120081 Bouchy et al. 120091 . It is fed by two op- 
tical fibers mounted at the focus of the 1.93-m telescope of the 
Haute-Provence Observatory (OHP, France). HAT-P-6 is bright 
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enough for observation in the high-resolution mode of the spec- 
trograph (A/AA = 75, 000) and with fast detector readout. 

The observations were carried out three days before full 
Moon. This was anticipated not to adversely impact the radial 
velocity accuracies, however, thanks to the large radial velocity 
shift between HAT-P-6 (-22.7 km s' 1 , N08) and the Moon (close 
to the barycentric Earth radial velocity, which in the direction of 
HAT-P-6 was 16.3 kms -1 during the observations). We used the 
second fiber input to measure the sky background, and confirmed 
that its radial velocity was always shifted by at least 30kms _1 
from that of HAT-P-6. We therefore did not subtract the sky 
background as this would have introduced additional noise. 

We collected 41 measurements during the transit night (JD 
2 455 430.5) between 20: 1 and 03 :40 UT. The sky was clear and 
the seeing stable around 1'.'6. The airmass first decreased from 
sec z = 1 .8 to 1 .0 then increased to 1 . 1 at the end of the night. We 
adjusted the exposure times between 500 and 800 s to maintain 
a constant signal-to-noise ratio (SNR) of 39 per pixel at 550 nm. 
Spectra of a thorium lamp were obtained at the beginning and 
end of the 7.5-hour sequence, as well as three times during the 
sequence. These wavelength calibration exposures show that the 
spectrograph drifted by ~ 2 ms _1 per hour. We linearly interpo- 
lated between the calibration points to correct for this drift. 

We used the SOPHIE pipeline (Bouchy et al. 120091 ) to extract 
the spectra, to cross-correlate them with a G2-type numerical 
mask, and to measure the radial velocities through Gaussian fits 
to the cross-correlation functions (CCFs) (Baranne et al. 1996; 
Pepe et al. 120021 ) . We discarded the four bluest orders of the 
SOPHIE spectra that have low SNR; this slightly improves the 
dispersion of the data around the model. Every spectrum pro- 
duces a clean CCF with a peak contrast of 1 1 .45 + 0. 10 % of the 
continuum and a full width at half maximum FWHM = 12.61 + 
0.08 kms _I . The photon-noise uncertainty of the radial veloc- 
ities, determined from the SNR of the spectra and the contrast 
and the FWHM of the CCF (see Boisse et al. l20T0l for the de- 
tails), is typically + 18 ms -1 . It is quite large by SOPHIE stan- 
dard because the star is early-type and a moderately fast rotator 
(Vsim, = 8.0 + 1.0 kms" 1 and 7/ eff = 6570 ± 80 K; see below). 

Figure [T] shows the SOPHIE velocities (see also the elec- 
tronic table) as well as the orbital solution and HIRES data from 
N08. Outside transit, the SOPHIE data are compatible (within 
the expected arbitrary offset) with the HIRES orbit. During tran- 
sit (lower panel), the SOPHIE data have obvious residuals rel- 
ative to a Keplerian orbit, with a pattern opposite that expected 
for an aligned, prograde transit (A - 0°), red-shifted radial ve- 
locity in the first half of the transit followed by a symmetric blue 
shift. Here, the observed pattern of blue-shifts during the first 
part of the transit and red-shifts in the second part, instead, is 
characteristic of the RM anomaly with a retrograde orbit. 

3. Analysis 

We fit our SOPHIE data together with the HIRES data from N08 
to determine the sky-projected obliquity in the HAT-P-6 system. 
Interestingly, two HIRES measurements were serendipitously 
obtained in transits, on December 16 th , 2006 and July 4 th , 2007 
(the two squares in the lower panel of Fig. [TJ. They agree with 
the SOPHIE measurements secured at the same orbital phases 
and thus with a retrograde orbit. 

The orbital elements of N08 determine the center of the 
August 21 st 2010 transit, T* = 2455 430.4563 + 0.0018 BJD. 
This +2.6-min uncertainty is negligible for our analysis. The 22 
SOPHIE measurements secured outside transit determine the 
systemic, barycentric stellar velocity as ysoPHiE = -22.594 + 
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Fig. 1. Radial velocity measurements of HAT-P-6 as a function 
of orbital phase. Filled diamonds represent the SOPHIE data 
(this paper) and open squares the HIRES data (N08). The lower 
panel zooms on the transit phases. The dashed line represents 
the Keplerian fit (ignoring the transit phases) while the solid 
line represents the final fit including the model of the Rossiter- 
McLaughlin anomaly. In each panel, the residuals of the final 
fit are plotted below the data. The vertical dotted lines show the 
times of mid-transit, first, second, third, and fourth contacts. The 
fitted values are V sin/ s = 7.5 + 1.6 kms -1 and A = 166° + 10°. 

0.007 km s _1 . The dispersion of their residuals from the orbital 
fit is 20.8 ms _1 . The error bars in the electronic table and Fig.Q] 
include an extra 9-ms~ 1 uncertainty added in quadrature to the 
photon noise to obtain a reduced x 2 of 1. Similarly the 13 
HIRES measurements outside transit determine a systemic ve- 
locity yniRES = -0.016 + 0.002 km s" 1 (HIRES velocities are 
relative) and have a 9.6-ms 1 dispersion around the Keplerian 
model. We quadratically added a 7.4-ms~ 1 jitter to their error 
bars (instead of the 8.6 m s -1 adopted by N08) to obtain ax 2 of 1 . 
The Can activity level measured from the SOPHIE spectra, 
log/?[j K = -5.03+0.10, agrees reasonably with log/?[j K = -4.81 
measured by Wright (120051 ). It predicts a lO-ms -1 stellar jitter 
(Santos et al.[2000 ) consistent with these added dispersion terms. 

We model the RM anomaly with the analytical first-moment 
approach developed by Ohta et al. ( 120051 ). We adopt the above 
values of the systemic velocities and mid-transit time, and the 
following parameters for the circular orbit taken from N08: ratio 
of the planetary and stellar radii R p /R t = 0.09338 ± 0.00053, 
scaled semi-major axis a/R* = 7.69 ± 0.22, semi-amplitude 
of the radial velocity variation K = 115.5 ± 3.5 ms~', or- 
bit inclination i = 85°51 + 0°35, and orbital period P = 
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3.852985 + 0.0 00005 days. We compute from a model atmo- 
sphere (Kurucz 1 19791 1 a linear limb-darkening coefficient e = 
0.57 + 0.10 in the 5300-6300 A wavelength range. 

The remaining two parameters of the Ohta et al. model are 
the sky-projected stellar rotational velocity V sin i s and the sky- 
projected spin-orbit angle A. We solve for these free param- 
eters through a grid search, and derive formal confidence in- 
tervals from the x 1 isocontours plotted in Fig. [2] (Hebrard et 
al. 120021) . The resulting values are Vsin/ S = 7.5 ± 1.2kms _1 
and A - 166°0 + 6°1. We separately propagate the effects of the 
uncertainties of the above fixed parameters on V sin i s and A, and 
find a + 1 . 1 km s" 1 effect on V sin i s , dominated by the uncertain- 
ties on a/R*, i, and ysoPHiE, and a ±8°5 effect on A, mainly from 
the ysoPHiE uncertainty. The other fixed parameters contribute 
negligibly to the A and V sin i s uncertainties. Quadratic addition 
of the two sources of uncertainties produces our final results: 
Vsinz' s = 7.5 + 1.6kms-' and A = 166° ± 10°. 
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Fig. 2. x 1 formal isocontours for our modeling of the Rossiter- 
McLaughlin effect as a function of A and V sin i s (see Sect. |3), 
using the formula from Ohta et al. (2005 ) and the system param- 
eters from N08. The diamond shows the best values. 



The resulting Vsim, agrees with the Vsin/j = 8.0 + 
1 .0 km s" 1 which we derive from the width of the CCF of HAT-P- 
6 (Boisse et al. 2010), as well as with V sin/, = 8.2 + l.Okms -1 
which N08 similarly derived from line broadening. The good 
agreement is somewhat surprising, since rotation broadens the 
lines of HAT-P-6 by significantly more than the SOPHIE reso- 
lution. As discussed by e.g. Hirano et al. (1201 Oat and Simpson 
et al. (120101 ). naive modeling of the RM anomaly could produce 
biased V sin i s measurements under such circumstances. 

Figure Q] plots the full dataset and the final fitted model. 
The SOPHIE and HIRES velocities have dispersions around this 
model of respectively 19.5 ms" 1 and 9.0ms" 1 . These are sim- 
ilar to the dispersions of the radial velocities measured outside 
transit around the Keplerian model, 20.8 m s" 1 and 9.6 m s" 1 (see 
above). Our RM model is thus a good description of the mea- 
surements, with a reduced ;f 2 of unity. 



4. Discussion 

The orbit of HAT-P-6b is retrograde with respect to the spin of 
its host star (Fig. [3), but the true, unprojected angle iff between 
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Fig. 3. Schematic view of the HAT-P-6 system from the Earth. 
The dotted and dashed-dotted lines represent the orbital angular 
momentum and the sky-projection of the stellar spin axis. The 
red and blue hues on the star surface indicate the Doppler shift 
from stellar rotation. The sky-projected obliquity A and its un- 
certainty are displayed in grey and green. 



the two angular momenta remains uncertain. We have only mea- 
sured the sky-projected spin-orbit angle, A = 166° + 10°, related 
to iff by cos if/ = sin i s cos A sin i„ + cos i s cos i (e.g. Fabry cky & 
Winn 2009 whose notations and definitions we adopt: i is the 
inclination of the planetary orbit and i s that of the stellar spin). 
The angle i Q is precisely known from the planetary transits but 
i s remains unknown. The star cannot be a priori assumed to 
be seen edge-on (i s i a 90°), especially here where A + 
demonstrates some misaligment. Comparisons of the measured 
V sin i s to estimate the stellar rotation period P rot and radius R* 
can in principle constrain i s , but here with ambiguous results. 

On one hand, the measured R' HK implies a short stellar ro- 
tation period, P rot = 3.5 days (Mamajek & Hillenbrand 2008). 
From/?* = 1.46 + 0.06 R G and VshU, = 8.0 + l.Okms" 1 , we 
obtain P 10t / sin i s = 9.2 + 1.5 days. Reconciling the two numbers 
requires an almost pole-on star, with i s =; 20° or ^Ofl implying 
iff 110°. Schlaufman (2010), on the other hand, predicts V ro t = 
9.1 km s" 1 from the age and mass of HAT-P-6. This is similar to 
the measured V sin i s , so would instead imply a more edge-on 
star, with i s 60° or 120°El and ip =* 145°. Neither approach 
is expected to be accurate. Here they give inconsistent estimates 
for P rot and therefore for i s and iff. Asteroseismology (Gizon & 
Solanki 2003), polarization of magnetic dipoles (Lopez-Ariste 
et al. 120 101 ). or accurate photometry may hopefully provide reli- 
able i s and P rot measurements. In the mean time the value of iff 
remains unknown, except for its being above 90°. 

HAT-P-6b is the seventh planet identified as having a retro- 
grade orbit, a fter W ASP-2b, 8b, 1 5b and 17b (Triaud et al. EoTOl 
Queloz et al. I20TU1 Bayliss et al. l20T0l ) a nd HA T-P-7b and 14b 
(Winn et al. |2009cl I2010al Narita et al. |2009l . Three planets 
additionally seem to have nearly polar orbits: CoRoT-lb (Pont 
et al. [20T0l ). WASP- lb (Sim pson et al. [20T0T >. and HAT-P-llb 
(Winn et al. 1201 Obi Hirano 1201 Obi ). The last of those is also 



1 The z, a* 200° and 340° formal solutions can be excluded since we 
know that the orbit is retrograde (\A\ > n — i a ). 

2 Here again, the retrograde orbit excludes i s ^ 240° and 300° . 
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likely retrograde. HAT-P-llb in addition is the only Neptune- 
mass planet with an obliquity measurement, with all other plan- 
etary systems with RM measurements being jovian mass planets. 

The processes that produce tilted systems remain a matter of 
debate (see Sect. [TJ. Whether a single mechanism explains all 
close-in planets or whether different processes produce aligned 
and misaligned systems is unknown. Even whether the planetary 
orbits acquired obliquity or the stellar spin instead acquired tilt 
is uncertain. Those mechanisms, in any case, cannot require too 
narrowly specified conditions, since tilted systems are common. 

Schlaufman (120101 1 and Winn et al. (2010c) hypothesized 
that misaligned planets preferentially orbit hot stars. The 6570 + 
80 K effective temperature of HAT-P-6 (N08) and the retrograde 
orbit of HAT-P-6b support this trend. As suggested by Johnson 
et al. (120091 1 and Hebrard et al. (2010), planetary mass appar- 
ently is also a key parameter. The seven retrograde and three 
polar planets all have masses below 3 Mj up . In that mass range, 
about a third of the planets with RM observations presents such 
extremes obliquities, with the other two thirds having A ^ 0°. 

Measured planets with masses above 4 Mj up are fewer but 
they seem to have a different obliquity distribution: four (over 
six) of these massive planets are misaligned, but none with ultra- 
high obliquities, i.e. none on retrograde orbit and even none with 
\A\ > 50°. Just two of the massive planets show A =* 0°, HAT-P- 
2b ( Winn e t al. l2007l Loeillet et al. l2008b and WASP- 18 (Triaud 
et al. 1201 01) . Curiously, the period distribution of exoplanets de- 
tected from radial velocity surveys changes at a similar charac- 
teristic mass. There is a distinct lack of planets with M p > 4 Mj up 
in short-period orbits (< 100 days), and radial velocity surveys 
detect no planets above that limit outside multiple star systems 

(Udry et al.|2002 Boisse et al.|20l3- 

The changes in the /l-distribution with planetary mass is il- 
lustrated in Fig.|4]by the 37 systems with published A and M p . 
The inset zooms on the Jupiter-mass planets, with their 2/3 of 
aligned and 1/3 strongly misaligned planets. The massive plan- 
ets, by contrast, mostly show significant but moderate misalign- 
ments. This suggests that distinct mechanisms dominate for dif- 
ferent planet masses, with a critical mass near 3.5 Mj up . If a com- 
bination of planet scattering, Kozai mechanism, and tidal circu- 
larization explains most oblique, close-in planets (Nagasawa et 
al. 12008 ; Morton & Johnson 2010), planetary mass is a natural 
control parameter. These conclusions however currently suffer 
from small-number statistics. Additional Rossiter-McLaughlin 
observations of massive planets are needed to confirm that most 
of them are tilted and none of them are polar nor retrograde. 
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